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Multifrequency radio continuum observations of 
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ABSTRACT 

We present high-sensitivity radio continuum observations with the VLA and Ryle 
Telescope at 1.5, 4.9, 8.4 and 15.4 GHz of the dwarf irregular galaxy NGC 1569. 
The radio data show an extended, irregularly-shaped halo with filamentary structure 
around the galaxy. The spectral index maps reveal an unusually patchy distribution 
with regions of flat spectral index extending into the halo. The data allow us to perform 
a spatially-resolved spectral fitting analysis of the continuum emission from which we 
derive maps of the thermal and synchrotron emission. The thermal radio emission is 
concentrated towards the brightest H ii region west of the super star clusters A and B 
whereas the distribution of the synchrotron emission peaks in a bar-like structure in 
the disk extending between the two clusters. The total flux density of the thermal radio 
emission allows us to derive the integra ted synchrotron spectrum and we confirm the 
break in the spectrum that was found bv llsrael fc de Bruvn I l)1988l) . We discuss various 
possibilities that could produce such a break and conclude that the only mechanism 
able to fit the radio data and being consistent with data at other wavelengths is a 
convective wind allowing cosmic ray electrons to escape from the halo. 

Key words: galaxies: individual: NGC 1569 - galaxies: starburst - radio continuum: 
galaxies - cosmic rays - convection 



1 INTRODUCTION 

NGC 1569 is a nearby d warf irregular galaxy (we adopt a 
distance of 2.2 ± 0.6 Mpc. Ilsraeill988ll . currently in the af- 
termath of a starburst. The high star-formation rate (SFR) 
has p roduced numerous young star-clusters (Hunter et al. 
I2OOOI) . the most lumino us being the super star clusters (SSC) 
A and B jAblejll97ih . The interstellar medium (ISM) re- 
flects the imprints of the active star formation in differ- 
ent ways. The atomic gas exhibits an unusually high ve- 
locity dispersion ( Stil & Israel 2000) and there is a re- 
ma rkable hole in the HI dis tribution, situated close to SSC 
A ilsrael fc van Drieill99(]H . The Ha emission is particu- 
larly unusual with a very filament ary structure, reminis- 
cent of an explosion IWalleil lTggl . and showing a com- 
plex vel ocity field consistent with th e presence of many 
bubbles jTomita. Ohta fc Saitdl 19941). From the kinemat- 
ics of the Hq filaments ijd ^Vaucouleurs et alj|l974 [Waller 



I991I. Heckman et al.ll995l . lMartin. Kobulnicky fc Heckmaii 



^002 ) and from X-ray emission ilHeckmar^^^T ' ^QSr 

Martin et alirioOSll the presence of an outflow of gas can 
be inferred. A possible trigger for the starburst is the inter- 
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action with a 7 x 1 0^ M© HI cloud at 5 kpc from NGC 1569 
JStn fc Israeill99dl . 

Numerous studies at different wavelengths have shed 
light on the star formation history. The analysis of colour- 
magnitude diagrams shows an increased level of star forma- 
tion (SF) starting about 100 Myr ago and stopping 4-10 
Myr ago (.Vallenari fc Bomanjri996llGreggio et al. igg&l . A 
drop of t he SFR about 4 - 10 Myr ago has also been in- 
ferred by iHeckman et alJ lll995ll from the paucity of stars 
more massive than about 20 - 25 Mq and from the low 
excitation temperature of the ionized gas derived from spec- 
troscopy. Assuming a Salpeter initial mass function (IMF), 
with a slope o f 2.35, and lower and upper mass cut-offs at 0.1 
and 120 Mrr. . ICreggio et al.l lll998D found a SFR of 0.5 Mq 
yr^^ with little change during the starburst phase (for an 
IMF slope of 2.6 the rate increases to 1 Mq yr ~^). This SFR 
is ver y high compared to ot her dwarf galaxies llGreggio et alJ 
Il998h . iHunter et all J2000h . on the other hand, found that 
most stellar clusters identified from HST observations had 
ages of less than 30 Myr, and suggested on this basis that 
most of the star formation took place towards the end of 
the starburst period . A similar conclusion was reached by 
I Anders et all (|2003) , who analysed multi- wavelength HST 
observations of a large number of star clusters in NGC 1569 
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and found from their derived age distribution that a major 
burst started about 25 Myr ago. 

Ahhough the SFR has dechned during the last 4-10 
Myr, it is still high, at least at some localized positions, as 
demonstrated by the large number of compact H ii regions 
JWaUeilll99l|^ and the young age of only a few Myr der ived 
for SSC A iGonzalez Dekado et alll997^ . IWalleI^ lll99lfl de- 
rives from the Ha data a present SFR of 0.4 M© yr~^ (as- 
suming an "e xtended Miller-S calo IMF" with a slope of 2.5 
above 1 Mq, |Kennicut3[l983). The inferred SFR would be 
reduced by a factor of 2 when taking into account only Ha 
emission from H ll regions (i.e. neglecting diffuse Ha emis- 
sion) and thus directly linked to massive star formation. This 
shows that the SFR has not decreased by more than a factor 
of about 3-5 during the last ~ 10 Myr compared to the last 
100 Myr. 

In this paper we present high-resolution radio contin- 
uum data at 4 wavelengths for NGC 1569. These data allow 
us to decompose the radio emission into its two components: 
free-free emission which is characterized by a spectral index^ 
of a = 0.1 and that originates principally in Hll regions, 
heated by stars with masses above about 20 M©, and syn- 
chrotron emission from diffuse relativistic cosmic ray elec- 
trons (CREs). The synchrotron spectrum is steeper, with 
a typically between 0.5 and 1.0, than that of the free-free 
emission and therefore dominates the radio spectrum at low 
frequencies, typically below 10 GHz. We assume that the 
main sites for the acceleration of cosmic rays (CRs) are su- 
pernova remnants (SNRs) associated with Type II and lb 
supernovae ( SNe). Such SNe have progenit ors with masses 
above 8 Mg iKennicutlill984 lBaronlll99 1 ) with life times 
of up to a few times lO'^ yr. After their injection into the 
ISM, CREs suffer energy losses mainly due to synchrotron 
and inverse Compton emission, which make their energy and 
radio spectra steepen in a characteristic fashion on typical 
time- scales around 10^ yr. 

Thus, the fitting of the radio continuum maps al- 
lows us to trace the sites of present SF (positions with 
free-free emission) and probe the time since the CREs 
were last accelerated, and hence previous epochs of SNe 
and SF. High frequency observations, preferably above 10 
GHz, are necessary for a reliable separation of the ther- 
mal and non-thermal radio emission. Previous studies with 
equivalent data sets were done for the starburst galaxy 
NGC 2146 (.Lisenfeld et al 199^) . the giant spiral galaxy 
NGC 1961 llL^enM d^^^|l99^) a nd the interacting gal ax- 
ies NGC 4490/4485 dClemens. Alexander fc GreeJl999ft . 



2 OBSERVATIONS AND DATA REDUCTION 

NGC 1569 was mapped with the VLA^ at 8.4 GHz and with 
the Cambridge Ryle Telescope (RT) at 4.9 and 15.4 GHz. 
Add itionally, existing VLA data (B-array) at 1.46 GHz 
JCondon 1985 ) were available. When the observations of 
NGC 1569 were made (between 1990 and 1993), the RT had 

^ In this paper we define the spectral index via S{i') oc 
where S{y) is the flux density at frequency v. 
^ The VLA is operated by the National Radio Astronomy Ob- 
servatory for Associated Universities Inc., under a cooperative 
agreement with the National Science Foundation. 



Table 1. Multifrequency radio maps 



Frequency Telescope r.m.s noise'^' 

[GHZ] [mjy/beam] 



1.5 


VLA (B-array) (2) 


0.12 


4.9 


Ryle 


0.13 


8.4 


VLA (B,C,D-array) 


0.09 


15.4 


Ryle 


0.15 



(1) at a resolution of 6 X 6 "(as the maps in Fig. 1) 

(2) from lCondoiil ^1983^ 



a band-width of 280 MHz split into 28 10-MHz frequency 
channels. Together with the minimum baseline of 18 m this 
provides excellent temperature sensitivity at both 4.9 and 
15.4 GHz. The maximum baseline of the RT is 4.6 km. 

The 1.5 GHz data were from ACondort .1983fl and no 
additional reduction was required. This map was made with 
uniform weighting and a Gaussian taper in the iM^plane of 
~ 30000 A, resulting in a ~ 6"synthesized beam. 

The 8.4-GHz VLA observations were in the B, C and 
D arrays and were observed during 1990/91. 3C286 and 
3C48 were used as flux calibrators. Reduction of the data 
followed standard VLA procedures in AlPS. The data were 
self-calibrated to convergence using a model consisting of 
the positive CLEAN components from the task CALIB. 

The RT observations were carried out between Septem- 
ber 1990 and February 1993. A total of five telescope con- 
figurations at each observing frequency were employed giv- 
ing a nearly fully filled aperture out to a baseline of about 
1 km. A minimum of two 12-hour runs in each configu- 
ration were obtained (any observation badly affected by 
interference was rejected). This gave a resolution of ap- 
proximately 1.0 X 1.0cosec((5) arcsec^ at 15.4 GHz and 
3.0 X 3.0cosec((5) arcsec^ at 4.9 GHz. At 4.9 GHz observa- 
tions of phase-calibrators were made at the beginning and 
end of each run, while at 15.4 GHz calibration observations 
were interleaved with those of NGC 1569. 3C 286 and 3C 48 
were observed regularly as flux calibrators. Calibration and 
data-editing (principally to remove narrow-band and time- 
varying interference) was performed in the MRAO package 
POSTMORTEM, with Subsequent reduction in AlPS and the 
MRAO package anmap. 

The maps at 4.9, 8.4 and 15.4 GHz were matched to 
the same spatial sensitivity as the 1.5 GHz map by restrict- 
ing the Mi>-range to values between lOOOA and 30000A, so 
that all maps were sensitive to spatial scales between about 
6" and 120". After this, the maps were CLEANed. AU maps 
have been corrected for the primary beam responses of the 
telescopes. This set of four images is shown in Fig. 1 and a 
summary of their characteristics is given in Table 1. 

The data at 1.5, 4.9 and 8.4 GHz are affected by the 
lack of short baselines which is reffected as reduced inte- 
grated ffux densities when compared to single-dish measure- 
ments (see Tab. 2). Only the map at 15.4 GHz, where most 
of the emission is very concentrated, does not show a con- 
siderable flux deficit. In order to make the spectral index 
comparison as accurate as possible, it is therefore crucial to 
ensure that the spectral index is determined from images all 
of which are sensitive as far as possible to the same range 
of spatial structure. The B-array VLA map at 1.5 GHz had 
the smallest range of spatial sensitivity and restricted there- 
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Figure 1. Maps of NGC 1569 at 1.5, 4.9, 8.4 and 15.4 GHz. The resolution is 6.0 X 6.0 arcsec'^ for each map. The contour levels start at 
two times the rms noise (see Table 1), corresponding to 0.24 mjy/beam, 0.26 mjy/beam, 0.18 mjy/beam and 0.30 mjy/beam at 1.5, 4.9, 
8.5 an d 15.4 GHz, respectively, and increase by factors of V2. In the 1.5 GHz map the positions of the strongest H II regions of lWalleii 
Jl99j) are shown, as well as the positions of the super star clusters (SSC) A and B. 



fore the uii-range used in the other maps. To achieve the 
match, the Whrange of the data in the aperture plane used 
in the construction of the images above 1.5 GHz was re- 
stricted during the mapping process so that the minimum 
and maximum basehnes (measured in wavelengths) used in 
the imaging were the same; this involved removing baselines 
from all the high-frequency data. While it is not possible 
to ensure that the precise sampling of the aperture plane is 
the same in all cases, for both the VLA and RT the dense 
sampling near the shortest baselines leads to similar overall 
coverage of the uv plane. 



emission lies beyond the optical galaxy in an irregularly- 
shaped halo. Prominent is an arm of emission to the West 
of t he galaxy wh ich coincides with a similar structure in 
Ha iWalleilll99 j) . Along the optical disk the radio emission 
shows three (at 15 GHz only 2) peaks. The two peaks vis- 
ible at all four wavele ngths coincid e with the two brightest 
Hll regions found bv lWalleil lll99lfl (number 2 and 7 in his 
nomenclature) . Apart from these features, the overall distri- 
bution of t he radio emis sion is in general very similar to the 
Ha map of lWalleil lll99lll . At higher frequencies the emission 
becomes more concentrated towards the major axis. 



3 THE RADIO PROPERTIES OF NGC 1569 

3.1 Structure of the radio emission 

Figure shows the radio emission at the four frequencies at 
a common resolution of 6 arcsec. At 1.5 GHz most of the 



3.2 Spectral index map 

The spectral index maps between 1.5 and 4.9 GHz, between 
1.5 and 8.4 GHz and between 8.4 and 15.4 GHz are shown 
in Fig. 13 Most striking is the extreme patchiness of these 
maps. The spectrum around the two main maxima of the 
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Table 2. Integrated flux densities 
Frequency Flux density Telescope 



Reference 



1 iVi.J.±Zj 


[mjy] 






38 


2300 ± 580 


CLFST 


1 


57.5 


1600 ± 500 


Clark Lake 


2 


151 


940 ± 180 


CLFST 


1 


610 


610 ± 20 


WSRT 


3 


1415 


440 ± 15 


WSRT 


3 


1415 


410 ± 25 


WSRT 


5 


1465 


285 ± 30 


VLA 


6,16 


1490 


411 


VLA 


7 


2695 


270 ± 50 


Effelsberg 


8 


2700 


365 ± 23 


Green Bank 


9 


4750 


263 ± 20 


Effelsberg 


10 


4850 


202 ± 19 


Green Bank 


11,12 


4919 


210 ± 20 


Ryle 


6 


4995 


280 ± 20 


WSRT 


3 


5000 


276 ± 42 


Green Bank 


9 


6630 


235 ± 30 


Algonquin 


13 


8415 


125 ± 12 


VLA 


6 


10700 


156 ± 16 


Effelsberg 


10 


10700 


154 ± 4 


Effelsberg 


14 


10700 


171 ± 10 


OVRO 


15 


15360 


116 ± 13 


Ryle 


6 


24500 


96 ± 8 


Effelsberg 


10 


The tables 


includes the 


data listed in 


Israel & 



(1988) (except for the points that 
able) as we ll as more recent data. 
11221), (2) Usraol & van Driel' (T 
198gli . (4) Icondon ct al.„. 1,199a) 
this w ork, (7)^^^^^^"^ 
\ 1989t). (9) ISulentid 



they rejected 
References; (1 



de Bruyn 
as unreli- 
jHowarth 



11 

r97e) 



3) [Israel fc ^e Bruvn 
I Jl98(t 



Hummel I 



(8) iPfleiderer. Durst fc Geblerl 
(10 ) iKlein fc Gravel il986l). 



(Ill iGregorv fc CondonI il99ih (12) iBecker. White fcEdwardsl 
J 199J). (13) iMcCucheonI il973^ . (14) iNiklas et al.1 Jl995l) . (15) 
llsrael van der Hulsj il98,'j) (16) ICondonl ^1983^ 
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Figure 2. Radio spectrum of NGC 1569. The data points are 
taken from Table 2. The triangles show our data and the crosses 
the data of other authors (shown as averages for observations at 
similar frequencies). The full line shows a single-power law fit 
with a = 0.47 to the data and the dashed line the thermal radio 
emission derived from the spectral fitting (sec Section 4). 




Figure 3. Spectral index distribution between 1.5 and 5 GHz 
(top), 1.5 and 8.4 GHz (middle) and 8.4 and 15.4 GHz (bottom) 
at a resolution of 6.0x6.0 arcsec^, overlaid over contours of the 
1.5 GHz emission (contour levels as in Fig. 1). The spectral index 
is calculated for all pixels where the flux exceeds three times the 
r.m.s. noise. The greyscales range from 0.0 to 0.8 (ai.s— sgHz and 
"8.4-i5GHz) and from 0.2 to 0.8 («i.5-8.4GHz)- 



radio emission and along the disk is flat and there is a trend 
for the low-frequency spectrum to steepen away from the 
galactic disk, up to values of a = 0.9 at the outermost end 
of the halo. This steepening does however not occur every- 
where and not in a regular way. There are numerous filamen- 
tary structures where the spectral index is much lower than 
in its surroundings and at some places the spectrum even 
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flattens with increasing distance from the centre. The high- 
frequency spectrum (a8.4-i5GHz shows an elongated region 
in the centre of the galaxy with a steep spectrum indicating 
the predominance of synchrotron emission in this region. 



4 SPECTRAL MODEL FITTING 

The radio-continuum emission consists of both a non- 
thermal (synchrotron) component and thermal emission. 
We use a spectral- fitting technique (th e same as in 
iLig gnfeld et alJ l|l99?j l. and similar to that in lClemens et all 
to estimate the two contributions and to determine 
spatial variations in the form of the synchrotron spectrum. 
For the set of maps of Fig. 1 we fit a function of the fol- 
lowing form at all pixels where the fiux densities at the four 
frequencies exceed three times the noise: 



I(u) = SoC{xB,i)+Tou- 



(1) 



where Tov~'^'^ is the form of the spectrum from an opti- 
cally thin Hll region, SoC{x-B,"i) is the synchrotron spec- 
trum from an aged electron population with injection index 
7 which have been radiating for a time t and = u/vb, 
where i/b is the break frequency (see below). During the fit- 
ting procedure the minimum in is found by varying So, vb. 
and To with the additional constraint that they are positive 
definite. 

The form of the synchrotron spectrum C(a;B, 7) depends 
on many processes, but to simplify the problem we make the 
following approximations: 

(i) The dominant electron energy loss mechanism is by 
synchrotron radiation and inverse Compton losses. These 
energy losses are proportional to the energy density of the 
magnetic field, C/b, and the energy density of the radiation 
field, ?7rad, respectively. The magnetic field can be estimated 
from the standard minimum energy assumption. With a ra- 
dio size of 3.5', a high-frequency cut-off of 8 GHz, a low- 
frequency cut-off of 100 MHz and a spectral index of 0.45, 
this yields B = 12/iG (Israel & de Bruyn 1988). C/rad can 
be estimated from the bolometric luminosity, Lboi, of the 
galaxy which we approximate here as the sum of the in - 
frared, optical and UV emission taken from 'Walleij il99ll) . 
(iboi = 7.4 X 10^^ W), by assuming that the radiation emit- 
ting matter is distributed spherically: 



c-kB? 



^ ^ X 1.4 X 10"^^ 



(i?/kpc 



eV 



(2) 



We derive, with R — 1.7 kpc (the average of major and 
minor axes of 3.6 and 1.8 arcmin), f7rad = 3.6 eV cm~"^. 

(ii) We neglect propagation effects within the galaxy. A 
spatially resolved model, as has been done for NGC 2146 
( Lisenfeld ct al. 1996), is not possible because of the extreme 
complexity of this galaxy making it impossible to derive the 
distribution of the sites of CR acceleration as a function 
of time. Given our small beam size (6 arcsec corresponds 
to only 66 pc) this assumption most likely does not repre- 
sent reality, because CREs are able to diffuse this distance 
(assuming a diffusion coefficient of lO'^^ cm^ s~^ as in our 
Galaxy) in about 10* yr, a time much shorter than the en- 
ergy loss time of typically lO'^ yr. However, this simplifying 
assumption still allows us to infer a realistic picture of the 



age distribution of the CREs, because the propagation in 
NGC 1569 is to a large extent driven by convection, as ar- 
gued in Sect. 5 and 6, and thus by an energy-independent 
process, which does not alter the shape of the synchrotron 
spectrum. Therefore, we expect that the CRE spectrum in 
each pixel mainly refiects the energy losses. 

(iii) We assume that the electrons which are radiating 
within a given resolution element were all accelerated (in- 
jected) a time t ago (the electron age). An alternative 
hypothesis for the s tar formation history was tested by 
IClemens et alj il999ft . They adopted in addition to the 
present star formation history a continuous injection of CRs 
starting some time t ago. The results for both star formation 
histories were practically identical. 

(iv) The electron loss time is much greater than the 
pitch angle randomization time — the so-called JP model 
( Jaffc & Pcrola 1974). 

(v) The electron injection index is 1.8, close to the value of 
2.0 predicted from simple theoretical models of shock accel- 
eration, and matching the observed low-frequency spectral 
index of NGC 1569. 

With these assumptions the electron energy distribution 
with an injection index 7 is given by: 



N{E,e,t) = 



NoE-'<{l - e{UB + U,^i)Ety~^ 




E <Eb 
e;>: Eb 



(3) 



with e = (4/3)crT/(7TioC), where ctt is the Thomson cross- 
section, and Eb = [^{Ub + t^rad)*]^^ is the break energy. 
The synchrotron spectrum is then obtained by convolving 
this electron energy distribution with the spectrum of the 
synch rotron emission of a single electron (see Lisciff cld et alJ 
il996r for explicit formulae). The break frequency is given 
by: 



UB = {3/2)ugEb oc 



B 



(t/rad + UB^t^ 



(4) 



where vq is the gyrofrequency. 



4.1 Results 

Maps of the synchrotron and thermal emissions and break 
frequency (which can be translated to CR electron age) are 
shown in Fig. |1| 

The integrated flux density of the thermal emission 
(Fig. m top) that we derive at 1 GHz is ~ 100 mjy. This 
number should be unaffected by the missing extended flux 
of our interferometer observations because the thermal emis- 
sion is very spatially concentrated as can be seen from the 
15.4 GHz map in Fig. The total thermal radio emission 
that we derive by our spectral fitting method corresponds 
very closely to the value of 97 mJy at 1 GHz estimated 
by [Israel & de Bruvn (1988) from reddening-corrected Ha 
measurements. Most of the thermal radio emission is con- 
centrated in the maximum of the radio continuum emission 
a nd coincides w ith the strongest Hll region (number 2) 
of lWalleJ 1I1991Y Correspondingly, the second strongest H 11 
region of Waller (number 7) is the second peak of the ther- 
mal radio emission. This good agreement between the Ha 
and thermal radio emission provides further support to the 
validity of our fitting method. 

The synchrotron emission (Fig. 2] middle) shows only 
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Figure 4. Greyscale maps of NGC 1569, showing the thermal 
emission (top), synchrotron emission (middle), both at 1 GHz, 
and break frequency (bottom) derived from the break in the radio 
spectrum, overlaid over contours of the 1.5 GHz emission. The 
greyscales are logarithmic in order to better highlight structures 
at different intensities. The range of break frequencies from 1 
(white) to 50 GHz (black) corresponds to electron ages from 1.4 X 
10® (black) to 1.0 X 10^ yr (white) for the adopted parameters of 
B and l/^ad- 



one pronounced peak at the south-eastern end, correspond- 
ing to the second strongest peak of the 1.5 GHz emission. 
Interestingly, this peak is not centered in the galaxy and is 
far away from the major sites of star formation traced by 
the thermal radio emission or the SSCs A and B. In the 
rest of the disk, the synchrotron emission is concentrated 
towards a bar-like structure along the major axis extending 
in between the positions of the SSC A and B. The peaks of 
the thermal and nonthermal radio emission do not coincide. 
At the maximum of the thermal emission there is only a 
moderate amount of synchrotron emission but surrounding 
this position there are several local maxima of synchrotron 
emission. 

The map of break frequency (Fig. 2] bottom) gives ages 
of several 10® years. The distribution is very patchy and can 
be divided approximately into two zones, one with relatively 
young CREs (ages of ~ 1 x 10® years) and other regions 
with older CREs (ages of 1 x 10^ years) . The young CREs 
coincide with the main synchrotron emission along the major 
axis of the galaxy, and extend slightly north of it. Towards 
the halo, the electrons generally become older although there 
remain patches of relatively young CREs. 

The large frequency range used - from 1.5 GHz where 
synchrotron emission contributes about 75% to the total 
emission up to 15.4 GHz where the thermal emission is the 
dominant (~ 90%) contribution - allows a reliable sepa- 
ration of both emissions based on their different spectral 
shapes. We therefore expect the error in the distribution of 
the synchrotron and thermal emission to be small and in the 
range of the errors of the individual maps. This is illustrated 
by a comparison of Fig. 3 and 4 which shows that the ther- 
mal emission is strong a places with a flat spectrum whereas 
strong synchrotron emission coincides with a steep spectrum 
up to high frequencies. The determination of the break fre- 
quency can be much less reliable at places with dominant 
thermal and little synchrotron emission and at those places 
where the synchrotron spectrum shows no clear curvature, 
leading to values of the break frequency below 1.5 or above 
15 GHz. Such values outside the range covered by our ra- 
dio data can only be estimates. The calculation of the elec- 
tron age from the break frequency depends directly on the 
adopted values for the magnetic and radiation field. Spa- 
tial variations of these two quantities, which are expected 
to be present but are not included in our model, produce 
therefore an error in the electron age. The separation of the 
thermal and synchrotron emission and the determination of 
the break frequency are however not affected by uncertain- 
ties in the magnetic and radiation field. 



5 TOTAL SYNCHROTRON SPECTRUM 

The thermal radio emission derived from the spectral fit- 
ting procedure of the radio maps allows us to determine the 
total thermal radio emission and, by subtraction, the total 
synchrotr on spectrum (Fig. El . Our r esult confirms the con- 
clusion of llsrael fc de Bruvn I (^8^ that the synchrotron 
spectrum shows an abrupt break, clearly visible in the data 
points above 10 GHz, because our derivation of the ther- 
m al radio emission co incides very well with the estimate 
of llsrael fc de Bruvn I dlQSSf) from extinction-corrected Ha 
emission. 
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Figure 5. The integrated synchrotron spectrum, obtained from 
the reliable data points in Fig.|^by subtracting of a thermal radio 
emission of Pthorm = 100 ^ (v /G}i7,)~^-^ mJy together with two 
fits to the data: time-dependent star-formation (dashed line) at 
t = 1 X 10^ yr (see Section 5.1 for details of the model) and 
convection into an open halo (full line, the model is described in 
Section 5.2). 

The shape of a synchrotron spectrum is determined by 
the energy losses, CR propagation processes and the star for- 
mation history. We have modeled various possibilities that 
could produce such a break and found that only two of them 
are in principle able to explain it. In the next subsections we 
will describe them as well as the other processes considered. 

5.1 Time-dependent star formation 

A temporal change in the SFR produces a change in the 
SN rate and therefore in the CRE production. If the time- 
scale of this change is shorter than the energy loss time scale 
of CREs it will b e noticeable in the synchrotron spectrum. 
iDeeg et alJ (^9^ have found evidence for this mechanism to 
be important in shaping the spectrum of several H ll galax- 
ies. NGC 1569 is a dwarf galaxy with variable SF in the 
recent past which can be expected to be reflected in its syn- 
chrotron spectrum. 

The break in the synchrotron spectrum is very sharp, 
therefore we will assume a change in the SFR as abrupt as 
possible. We assume that SF started at t = and ended at 
tend- The exact value of tend is not important, as long as 
it is sufficiently large (here we chose t^nd = 10* yr) to have 
achieved a steady-state situation for massive stars and CREs 
by that time. We further assume that the stars are produced 
according to a Salpeter-IMF and that all stars above 8 M0 
end their lives in a SN explosion which accelerates CRs. As 
in the spectral fitting, the main energy losses are assumed to 
be synchrotron and inverse Compton and we use the values 
for the magnetic field B and the radiation energy density 
derived in the previous section. 

The synchrotron emission is then calculated by solving 
the time-dependent equation for the electron particle density 
f(t,E): 



This equation takes into account CR acceleration in SN 
remnants with a source spectrum as (E / mc?)~'^ qsnvsT<i{t) 
(where vs^(t) is the supernova rate, calculated from the 
above defined SFR(t), and gsN is the number of relativis- 
tic electrons produced per supernova and per unit energy 
interval), and the subsequent radiative energy losses of the 
CR electrons (6 oc {Ura.d +U-b))- 

The solution of Eq. (5) is given by: 

/(£,*)= (^^^ qsN dt'iysNit')(^l-Eb{t-t')^ ,(6) 

where to = max{0,t — rioss) and Tioas(-E,t) — (Eb)^^ oc 
((C^rad + U-b)E)~^ is the life-time of the CREs against radia- 
tive energy losses. Finally the total synchrotron spectrum is 
obtained by convolving f{E, t) with the synchrotron emis- 
sion spectrum of a single electron. 

We calculated the spectra as a function of time and 
selected the best-fit spectra. There are two possible times 
when a break in the spectrum is produced: 

1) A good fit to the data can be achieved soon after 
the beginning of the starburst, at f = 1 x 10^ yr (see Figl^i, 
with a realistic injection spectral index 7 = 2.0. The break 
in the spectrum divides the low frequency range, in which 
the electrons have not suffered energy losses yet, from the 
high frequencies where the spectrum is steepened by 0.5 as a 
result of energy losses. The position of the break is given by 
the life-time of the electrons against radiative losses, which 
is 5 X IQp yr at 5 GHz (roughly the position of the break) 
with o ur values for the magnetic field and f/rad . 

2) llsrael fc de Bruvn I lll988l) suggested that the break 
in the synchrotron spectrum was due to the abrupt end of 
the CR injection about 5 x 10® yr ago. A cut-off at high ener- 
gies is then produced above which electrons have already lost 
their energy since the time of the last CR injection. We have 
tried to fit the data for times after the end of the starburst. 
Although a break is indeed produced by this mechanism, we 
could not find a satisfactory fit to the data. Two reasons 
are responsible for this failure: 1) A visible break cannot 
be produced at the high frequencies where it is observed. 
The main reason is that, even for an abruptly ending SFR, 
the SN rate does not cease immediately. If stars above 8 
M0, with life-times of 3 x 10^ yr, become SN, then it takes 
about this time until the SN rate has ceased after the end 
of the SF. This time-scale is longer than the CR loss-time 
scale of ~ 5 X 10® years of CR electrons emitting in the 
range around 5 GHz where the break is observed, so that 
these high-energy CREs do not experience an abrupt end 
in the injection. If the SN rate ended abruptly, the break 
in the spectrum would be more pronounced in better agree- 
ment with the data. 2) The low-frequency spectrum is very 
flat, requiring an unrealistically low injection spectral index 
(7 = 1.0) to explain it. 

5.2 Escape from the halo 

An alternative mechanism causing a break in the spectrum 
is escape of CREs from the galaxy. If the escape time-scale 
is comparable to the energy loss time-scale rioss a break of 
Aq = 0.5 is produced in the spectrum which is sufficient to 
explain the data if the break is sharp enough. Convection as 
the mode of propagation is required for this. If diffusion is 
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the dominant propagation mechanism, the frequency varia- 
tion of the expected radio spectral index is too smooth to 
account for the sharp break in the data (see next section). 

We model this situation as follows: CREs that are pro- 
duced in the main disk of NGC 1569 convect into the halo 
with a constant convection velocity, Vc, experiencing syn- 
chrotron and inverse Compton losses at the same time. 
When they reach a distance Zhaio from the disk, they can 
escape freely. The one-dimensional convection-energy loss 
equation describing this situation is: 

with 8{z) being the one-dimensional 5-function. The solution 
of this equation is given by: 

The synchrotron emission is calculated by first integrat- 
ing this emission from z = Q to zOmm(2max, Zhaio), where 
2max = Vc/{bE) IS the maximum distance that CREs can 
travel before having lost their energy to below E, and then 
convolving the result with the synchrotron emission spec- 
trum of a single electron. 

We use the above values of C/rad and Ub, and estimate 
-Zhaio = 0.7 kpc from the extent of the radio emission at 1.5 
GHz. Then, we obtain a good fit to the data with Vc ~ 150 
km s"^ (Fig.EJ. 

5.3 Other possibilities 

We have also tested other mechanisms that produce breaks 
in the spectral index, but none is able to produce the sharp 
break seen in NGC 1569. 

Escape from the halo, as described in the last subsec- 
tion, in combination with diffusion as the dominant prop- 
agation mechanism, produces too weak a frequency varia- 
tion of the radio spectral index. The reason for this dif- 
ference is that the scale length for diffusion (i.e. the dis- 
tance over which CR electrons can propagate before losing 
their energy), Xdiff{E) cx {D{E) /\Eh)Y^'^ , depends only 
weakly on energy, even for a constant diffusion coefficient 
D{E) = Dq. If the diffusion coefficient depends on energy 
as D(E) = -Dq-E^' with fi ~ 0.5, as derived for the Galaxy 
jBerezinskv et alJll990h . the energy dependence of A is re- 
duced even more. In the case of convection the scale length 
has a stronger energy dependence, XconviE) oc {V/{Eb)), 
and therefore produces a sharper break. 

In a closed halo model, a break in the spectrum can 
be caused by different energy loss processes that are impor- 
tant for different ranges of CRE energies. A break in the 
synchrotron spectrum of Aq — 0.5 can be caused by adia- 
batic and bremsstrahlung losses in addition to synchrotron 
and inverse Compton losses. In some spiral galaxies this can 
explain the flatten ing of the synchrot ron spectrum towards 
lower frequencies JPohl et al. lll991aD . Taking additionally 
into account ionization losses, that are important at low fre- 
quencies, a total break of Aa = 1.0 can be produced. This 
process is able to explai n the low-frequency flattening of the 
spectrum of NGC 4631 jPohl et alJll991bl) . We have tested 
these possibilities and even in the latter case, the spectrum 



has much too weak a curvature and is therefore not able 
to explain the sharp bend in the synchrotron spectrum of 
NGC 1569. 

A further possibility to explain a flat low-frequency 
spectrum are absorption processes affecting the synchrotron 
emiss ion. Whereas syn chrotron self-absorption is not impor- 
tant iPeeg et alJll993^ . free- free absorption can play a role. 
We have tested both absorption by a foreground layer and 
absorption of ma t erial mixed with th e emitting material (see 
iDeeg et al.lll993l . iPohl et al.l|l991bft and found that in nei- 
ther case this process is able to explain the sharp break in 
the spectrum. 



6 DISCUSSION 

6.1 Star formation history 

The maps of thermal radio and synchrotron emission trace 
the star formation history on different time-scales: Whereas 
the thermal radio emission traces the star formation of the 
last ~ 10^ years, the synchrotron emission probes a longer 
timescale. The lifetime of CREs in NGC 1569 range, with the 
parameters for B and Kad adopted here, between a few up to 
about 10 Myr for electrons radiating in the GHz range. This 
gives, added to the lifetime of SN progenitor stars of up to 
30 Myr, a span of 40 Myr that is probed by the synchrotron 
emission. 

The image of the thermal radio emission in Fig.|31shows 
that the recent SF is concentrated to a large extent in 
one peak, coinciding with the Hll region 2. This position 
furthermore coincides with the peak of the dust emission 
(Lisenfeld et al. 2002) an d is about 10 arcse c north of the 
peak of the CO emission llTavlor et al.lll999^ showing that 
this whole area is active in star formation. 

The peak of the synchrotron emission in the south-east 
shows a region that has been the site of star formation in 
the past. We can derive the SFR from the integr ated flux 
of this region (8 mjy) and derive, using eq. 21 in ICondon I 
tl992J, 8.3 X 10"'' Mq/jv for stars more massive than 5 Mq. 
This is a factor 24 less than the SFR derived from the peak 
of the thermal radio emission (20 mJy and 2 x 10"^ Mo/yr, 
using eq. 23 in lCondon I il992h ') showing that the SFR in the 
synchrotron peak has been less important than the present- 
day SF in the west of the galactic disk. Within the disk, the 
SF traced by the synchrotron emission has a bar-like shape 
and extends in between the SSC A and B. The ages of these 
clusters - the stellar populations are estimat ed to have ages 
up to 9 Myr iGonzalez Delgado et al.lll997l) - are younger 
than the SF traced by the synchrotron emission. This might 
indicate that the formation of the SSC has been triggered 
by the preceding epoch of close-by SF. 

The break in the integrated synchrotron emission can 
be explained by the effect of variations of the SF in the past. 
We have shown that it cannot be due the end of the star - 
burst 5 Myr ago, as suggested bv llsrael fc de Bruvn I Jl98d) . 
The break could only be explained by the beginning of the 
starburst about 10 Myr ago. This explanation is however not 
in agreement with several other studies that agree that the 
SFR has declined during the last 5-10 Myr and that the star- 
burst preceding this decline had a duration of abo ut 100 Myr 
JVallenari fc Bomansl[l99l ICreggio et alJll998ll . Therefore 



Convective wind in NGC 1569 9 



we conclude that the abrupt beginning of the SF some 10 
Myr ago, necessary to explain the break in the synchrotron 
spectrum, is unlikely to have taken place. 

6.2 Outflow from the halo 

Since varying SF has been excluded as a cause of the break 
in the synchrotron spectrum, the only possible explanation 
is escape of the CREs from the halo. Such an outflow is 
not unexpected: massive star formation and the related en- 
ergy input by SN explosions can lead to an outflow of mat- 
ter from the galactic disks into the halo. In starbursting 
dwarf galaxies like NGC 1569 this process is expected to be 
particularly important due to their relatively small gravita- 
tional potential. The unusually large size of the radio halo 
of N GC 1569 is by itself a sign of a h igh energy input by 
SNe JPahlem. Lisenfeld fc GoUall 1995t) . Further strong evi- 
dence for the presence of a galactic wind is the weak steep- 
ening of the spectral index in the halo indicating a con- 
vective outflow of CRs along vertical magnetic fleld lines 
which enables a fast (with respect to the energy loss time 
scale of the CREs) transport of CRs into the halo. Such a 
magnetic fleld perpendicular to the disk can be produced 
by stretching the fleld along the flow direction of a galac- 
tic wind (Lcsch ct al. 1989 ). It has been found in the star- 
burst galaxy M 82 llReuter et alj[l992f) . Recent polarized 
radio continuum observations at 8.5 GHz (VLA) and 10.4 
GHz (Effelsberg) show i ndeed an out wardly directly mag- 
netic fleld in NGC 1569 jMuhlell20o3l . supporting the pos- 
sibility of a convective outflow from NGC 1569. 

Further support for the outflow hypothesis is a general 
similarity between the radio continuum map at 1.5 GHz and 
the Ha image, where the irregular halo to the south and 
a long arm extending from the northwest of the galaxy is 
clearly visible in both images. At this low frequency the ra- 
dio emission is predominantly non-thermal. An overall sim- 
ilarity between the Ha and r adio emission in the halos has 
also been found in NGC 891 dPahlem. Dettmar fc HummeJ 
I 1994h . NGC 4631 and NGC 5775 l)Dettmailll992h . Both 
emissions are related to regions of active st ar formation 
within the galactic disk jPahlem et alJ^1994^ . A plausible 
explanation for this coincidence is an outflow of ionized gas 
due to the high star formation activity which at the same 
time drags the magnetic fleld along and enables therefore a 
convective outflow of the CREs ( I )alil(Mii d al. _li)95 , 2001). 

From our flt to the total synchrotron spectrum we de- 
rived a convection velocity of 150 km s~^. This value is 
higher t han the escap e velocity of this galaxy of 80 - 110 
km s~^ (lMartinlll993 ^ suggesting that the outflow is able 
to eject gas from the galaxy. An alternative estimate for 
the convective velocity can be made from the spectral index 
maps (Fig.laJ. In some fllamentary structures the spectral 
index does not steepen outwardly, indicating that the CREs 
do not suffer substantial energy losses. Assuming that the 
CR electrons were produced in the disk and taking an aver- 
age length of such fllaments in the maps of ai.5_5GHz (where 
the contribution by the thermal radio emission is still small) 
of 20 - 30" (corresponding to 220 - 330 pc) we can derive 
from the energy loss time-scale of 4.5 x 10^ yr at 5 GHz a 
lower limit for the convective velocity of 50 - 70 km s~^, 
consistent with the value derived from the total spectrum. 

Observations of the hot gas in the halo of NGC 1569 



provide independent evidence. iHeckman et alJ ^1995^ have 
found indications for a supernova-driven outflow from 
ROSAT X-ray and Ha emission. They flnd that about 25% 
of the Ha is in a large (~ 2 kpc) complex system with ra- 
dial velocities of as much as 200 km s~^. This velocity agrees 
well with our convection velocity derived from the total spec- 
tral index. A comparison of the Ha shell speeds to the es- 
cape velocity estimated from the rotation curve and taking 
into account a dark matter halo, Martin (1998) showed that 
NGC 1569 is likely to present a substantial loss of shell ma- 
terial from the galaxy. The velocity of the hot gas traced 
by X -r ay emission is well ab ove the escape velocity ijMartinI 
ll999^ . lMartin et alj d 20021) confirmed this conclusion with 
deep Chandra spectra and found that the outfiow is strongly 
metal-enriched, carrying nearly all the metals ejected by the 
starburst. 



7 CONCLUSIONS AND SUMMARY 

We have present high-sensitivity radio continuum observa- 
tions at 1.5, 4.9, 8.4 and 15.4 GHz of the dwarf irregular 
galaxy NGC 1569 at a high resolution (6 arcsec). These data 
allow us to perform a spatially-resolved spectral fitting anal- 
ysis of the continuum emission from which we derive maps of 
the thermal and synchrotron emission. The total flux den- 
sity of the thermal radio emission allows us to derive the 
integrated synchrotron spectrum. The results from our data 
and its analysis are: 

• The radio data show an extended, irregularly shaped 
halo with fllamentary structures around the galaxy. The 
spectral indices between 1.5 and 4.9 GHz and between 1.5 
and 8.4 GHz show an unusually patchy distribution with 
regions of flat spectrum extending into the halo. 

• The distribution of the thermal and synchrotron emis- 
sion trace the star formation history on different time-scale 
(~ lO'^ yr for the thermal emission and ~ 4 x 10^ yr for 
the synchrotron emission) in NGC 1569. The thermal radio 
emission is mainly concentrated in the brightest H ll region 
west of the super star cluster A and B whereas the syn- 
chrotron emission peaks in a bar-like structure in the disk 
extending in between the two star clusters. The super star 
clusters are very young (< lO'' vr. IConzalez Delgado et alJ 
'199/) and formed after the SF period traced by the syn- 
chrotron emission. It is thus possible that their formation 
has been triggered by the older SF in the bar. 

• We conflrm the break in the synchrotron spectrum that 
was found by Israel & dc Bruvn (1988). We discuss vari- 
ous possibilities that could produce such a break and flnd 
only two that are consistent with the radio data: an abrupt 
beginning of the star formation about 10^ years ago or a 
convective wind transporting CREs into an open halo. An 
abrupt end of a sta r burst several Myr ago, as suggested by 
llsrael fc de Bruvn I l)l988|) . is not able to explain the data 
because even in this extreme case the SN rate would not 
have declined immediately but only over a period of about 
3 X 10^ years. This time scale is longer than the energy loss 
time scale of CREs so that a sharp break would be consid- 
erably smoothed. An abrupt beginning of the starburst lO'^ 
years ago is not consistent with data at other wavelengths 
from which a decline (and not a n increase) in the SFR over 
this time scale is derived (e.g. IVallenari fc BomansI Il99^ . 



10 U. Lisenfeld et al. 



iGreggio et alJlQQ^) . We thus conclude that the only mecha- 
nism able to fit the radio data and being consistent with data 
at other wavelengths is a convective wind allowing CREs to 
escape from the halo. This is in agreement wi th X-ray studies 
indicating the presence o f a galactic wind jHeckman et alJ 
llQQfll.lMartin et al.ll 20o3^ . 



ACKNOWLEDGMENTS 

We would like to thank the referee, U. Klein, for useful sug- 
gestions which helped us to improve the paper. UL acknowl- 
edges support by the Spanish MCyT AYA2002-03338 and 
the Junta de Andaluci'a (Spain). 



REFERENCES 

Abies, H. D., 1971, Publ. U.S. Naval Obs. Sec. Ser. XX 
(IV), 61 

Anders, P., de Grijs, R., Fritze-v. Alvenslegen, U., Bissantz, 

N., 2004, MNRAS, 347, 17 
Arp, H. C, Sandage, A. R., 1985, AJ, 90, 1163 
Baron, E., 1992, MNRAS, 255, 267 

Berezinsky, V.S., Bulanov, S.V., Ginzburg, V.L., Dogiel, 
V.A., Ptsuskin, V.S., 1990, Astrophysics of Cosmic Rays, 
North Holland:Amsterdam 

Becker, R. H., White, R. L., Edwards, A. L., 1991, ApJS, 
75, 1 

Clemens, M.S., Alexander, P., Green, D.A., 1999, MNRAS, 
307, 481 

Condon, J. J., 1983, ApJS, 53, 459 
Condon, J. J., 1987, ApJS, 65, 485 

Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. F., 

Perley, R. A., Taylor, G. B., 1998, AJ, 115, 1693 
Condon, J. J., 1992, ARAA, 30, 575 

Dahlem, M., Dettmar, R.-J., Hummel, E., 1994, A&A, 290, 
384 

Dahlem, M., Lisenfeld, U., GoUa, G., 1995, ApJ, 444, 119 

Dahlem, M., Lazendic, J. S., Haynes, R. F., Ehle, M., Lisen- 
feld, U., 2001, A&A, 374, 42 

Dettmar, R.J., 1992, Fund, of Cosmic Phys., 15, 143 

Deeg, H.J., Brinks, E., Duric, N., Klein, U., Skillman, E., 
1993, ApJ, 410, 626 

de Vaucouleurs, G., de Vaucouleurs, A., & Pence, W., 1974, 
ApJ, 194, L119 

Gonzalez Delgado, R.M., Leitherer, C., Heckman, T., 
Cerviiio, M., 1997, ApJ, 483, 705 

Greggio, L., Tosi, M., Clampin, M., de Marchi, G., Lei- 
therer, C., Nota, A., Sirianni, M., 1998, ApJ, 504, 725 

Gregory, P. C, Condon, J. J., 1991, ApJS, 75, 1011 

Heckman, T.M., Dahlem, M., Lehnert, M.D., Fabbiano, G., 
Gilmore, D., Waller, W.H., 1995, ApJ, 448, 98 

Howarth, N. A., 1990, PhD thesis. University of Cambridge 

Hummel, E., 1980, A&AS, 41, 151 

Hunter, D.A., O'ConneU, R.W., Gallagher, J.S., Smecker- 

Hane, T.A., 2000, AJ, 120, 2383 
Israel, F. P., 1988, A&A, 194, 24 
Israel, F. P., van der Hulst, J. M., 1983, AJ, 88, 1736 
Israel, F. P., de Bruyn, A. G., 1988, A&A, 198, 109 
Israel, F. P., Mahoney, M.J., 1990, ApJ, 352, 30 
Israel, F. P., van Driel, W., 1990, A&A, 236, 323 



Jaffe, W.J., Perola, G.C 1974, A&A, 26, 423 

Kennicutt, R.C, 1983, ApJ, 272, 54 

Kennicutt, R.C, 1984, ApJ, 277, 361 

Klein, U., Grave, A., 1986, A&A, 161, 155 

Lesch, H., Crusius, A., Schlickeiser, R., Wielebinski, R., 

1989, A&A, 217, 99 
Lisenfeld, U., Alexander, P., Pooley, G.G., Wilding, T., 

1996, MNRAS, 281, 301 
Lisenfeld, U., Alexander, P., Pooley, G.G., Wilding, T., 

1998, MNRAS, 300, 30 
Lisenfeld, U., Israel, F. P., Stil, J., Sievers, A., 2002, A&A, 

382, 860 

Martin, C.L., 1998, ApJ, 506, 222 
Martin, C.L., 1999, ApJ, 513, 156 

Martin, C. L., Kobulnicky, H. A., Heckman, T. M., 2002, 

ApJ, 574, 663 
McCucheon, W.H., 1973, AJ, 78, 18 
Miihle, S., 2003, PhD thesis. University of Bonn 
Niklas, S., Klein, U., Braine, J., Wielebinski, R., 1995, 

A&AS, 114, 21 
Pfleiderer, J., Durst, C, Gebler, K.H., 1980, MNRAS, 192, 

635 

Pohl, M., Schlickeiser, R., Hummel, E., 1991a, A&A, 250, 
302 

Pohl, M., Schlickeiser, R., Lesch, H., 1991b, A&A, 252, 493 
Renter, H.-P., Klein, U., Lesch, H., Wielebinski, R., Kron- 

berg, P.P., 1992, A&A,250, 302 
Stil, J., Israel, F.P., 1998, A&A, 337, 64 
Stil, J., Israel, F.P., 2000, A&A,392, 473 
Sulentic, J. W., 1976, ApJS, 32, 171 
Tomita, A., Ohta,K., Saito, M., 1994, PASJ, 46, 335 
Taylor, C.L., Hiittemeister, S., Klein, U., Greve, A., 1999, 

A&A, 349, 424 
Vallenari, A., Bomans, D. J., 1996, A&A, 313, 713 
Waller, W. H., 1991, ApJ, 370, 144 



